Synthesizing the Lii attractor by parameter-switching 
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Abstract 



j> In this letter we synthesize numerically the Lii attractor starting from the generalized 

Lorenz and Chen systems, by switching the control parameter inside a chosen finite set 

of values on every successive adjacent finite time intervals. A numerical method with 

O fixed step size for ODEs is used to integrate the underlying initial value problem. As 

t-h numerically and computationally proved in this work, the utilized attractors synthesis 

. !_h algorithm introduced by the present author before, allows to synthesize the Lii attractor 

><! 

5— i starting from any finite set of parameter values. 
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1 Introduction 



Consider the following unified chaotic system (bridge between the Lorenz and Chen systems) 



Lii et aL 2002 



xi = (25p + 10)(x 2 -xi), 

x 2 = (28 - 35p) Xl + (29p - l)x 2 - Xl x 3 , (1) 
x 3 = xix 2 — (8 +p) /3x 3 , 

where the parameter p E [0, 1] . As it is known now, for p E [0, 0.8) ([I]) models the canon- 
ical Lorenz system |Chelikovsky fc Chen, 2002 , for p = 0.8 the system becomes Lii system 
Lii fc Chen, 2002] , while whenp E (0.8, 1], the system becomes Chen system Chen Ueta, 1 999 



Therefore, this system is likely to be the simplest chaotic system bridging the gap between 
the Lorenz and the Chen systems. 

The above three systems share some common properties such as: they all have the same 
symmetry, dissipativity, stability of equilibria, similar bifurcations and topological structures 
and belong to the generalized Lorenz canonical family |Chelikovsky fc Chen, 20 02 1. 

In the mentioned references, a positive answer to the question as if it is possible to realize 
a continuous transition from one to another system is given. 

In this letter, we present a discontinuous transition algorithm between the Lorenz and 
the Chen systems with whitch the Lii attractor can be synthesized. For this purpose, the 
parameter switching method introduced in |Danca et aL, 2008] is utilized. 

The present work is organized as follows: Section 2 presents the synthesis algorithm, while 
in Section 3 the Lii attractor is synthesized in both deterministic and random ways via the 
mentioned synthesis algorithm. 
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2 Attractors synthesis algorithm 

Consider a class of dissipative autonomous dynamical systems modeled by the following initial 
value problem: 

S : x = f p (x), x(0) = x , (2) 
where peK and f p : lR n — > M. n has the expression 

f p (x) =g(x)+ P Mx, (3) 

with g : M. n — > M. n being a vector continuous nonlinear function, M a real constant n x n 
matrix, xq G W 1 , and the maximal existence interval / = [0, oo). 
For the Lii system ([lp, one has: 

/ 25 ^ 

M= -35 29 , g(x) = (10 (x 2 - x\) , x 2 - xix 3 + 28xi,xix 2 - 8/3x 3 ) T 
\ -1/3 y 
with divergence div f p (x) < for p G [0, 1], so the system is dissipative. 

The existence and uniqueness of solutions on the maximal existence interval / are as- 
sumed. Also, without any restriction, it is supposed that corresponding to different p, there 
are different global attractors. Because of numerical characteristics of the attractors synthesis 
(AS) algorithm and for sake of simplicity, by a (global) attractor in this letter one under- 
stand without a significant loss of generality, only the approximation of the u-limit set, as 
Foias & Jolly, 1995 , is plotted after neglecting a sufficiently long period of transients (for 



in 



background about attractors see |Milnor, T 985|). 

Notation 1 Let V CR be the set of all admissible values for p and A the set of all corre- 
sponding global attractors, which includes attractive stable fixed points, limit cycles and chaotic 
attractors. Also, denote by Vn a finite subset of V for some positive integer N > 1 and the 
corresponding subset of attractors An C A. 



Because of the assumed dissipativity, A is a non-empty set. Therefore, following the above 
assumptions, a bijection between V and A, F : V — > A, can be considered. Thus, to each 
p G V corresponds a unique global attractor A p e A and conversely for each global attractor 
there exists a unique parameter value p EV. 



In Danca et al., 2008 , it is proved numerically that switching, indefinitely in some pe- 
riodic waythe parameter p inside Vn over finite time subintervals, while Q is integrated 
with some numerical method for ODEs with fixed step size h, any attractor of An can be 
synthesized. For a chosen N, consider / being partitioned in to consecutive sets of iV finite 
adjacent time subintervals Ij, : I — (I\ U I2 U . . . U In) U (1% U I2 U . . . U In) U ... of lengths 
Ati, i = 1, 2 . . . , N. If, in each subinterval Ii, while some numerical method with single fixed 
step size h integrates (J2|, p is switched as follows: p = pi, for t e Jj, Then, a synthesized 
attractor, denoted by A*, can be generated. The simplest way to implement numerically the 
AS algorithm is to choose Ati as a multiple of h. Thus, the AS algorithm can be symbolically 
written for a fixed step size h as follows: 

[mipi, TO2P2, • • • , ttinPn], (4) 

where are some positive integers (weights) and by m^pk one understands that in the k-th 
time subinterval Ik, of length rrikh, p receives the value pk- 

In |Danca et al, 2008] , it is proved numerically that A* is identical^ to A p * for 



JV 

E m kPk 

p* = ^ ■ ( 5 ) 

fc=i 



^denticity is understood in a geometrical sense: two attractors are considered to be (almost) identical if 
their trajectories in the phase space coincide. The word almost corresponds to the case of chaotic attractors, 
where identity may appear only after infinite time. Supplementarily, Poincare sections and Haussdorf distance 
between trajectories are utilized to underline this identity. 



For example, the sequence [lpi,2p2] means that vii\ = 1, m 2 = 2 and the synthesized 
attractor A* is synthesized as follows: in the first time interval I\ of length Ati = h, the 
numerical method solves ^ with p = p%] next, for the second time interval I2 of length 
At 2 = 2/i, p = P2, and the algorithm repeats. If we apply this scheme to ([T]) for p\ = 0.8 
(chaotic Lii attractor) and P2 = 0.959 (chaotic Chen attractor), one obtains the synthesized 
regular Chen attractor A* which is identical to A p * with p* = (p 1 + 2p 2 ) /3 = 0.906 in ([5]), 
corresponding to a stable periodic limit cycle. In Fig. 1, to underline the identity, phase plots, 
time series, histograms and Poincare sections superimposed were utilized beside Haussdorf 
distance between the two attractors which is of order 10 -2 -j- 10 -3 conferring a good accuracy 
to AS. 

It is noted that the AS algorithm can be applied even in some random way: because 

N N 

p* in (5) is defined in a convex manner (if denoting = m^j < 1, then p* = otkPk-, 

k=l k=l 

N 

with otk — 1) and based on the bijective function F, any synthesized attractor A* is located 

k=l 

inside the set An (all elements, i.e. attractors, are ordered with the order endowed by F) and 



whatever (random) scheme ^ is used, the result is the same Danca, 2008 



The random AS can be implemented e.g. by generating a sequence §4§ with a random uniform 
distribution of p |Danca, 2008| which is supposed to generate all the integers 1, . . . , N (Fig. 
2). 

Now, p* is given by the following formula: 

N 

Y m 'iPi 

P* = ^ (6) 

1=1 

where m! i counts the number of pi . Obviously, now, I has to be chosen large enough, such 
that ^ can converge to p* (the precise value in this case for p* could be obtained only for 
I = [0, 00) ). 
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Remark 2 i) The AS algorithm is useful in the applications where some p are not directly 
accessible. 

ii) The AS algorithm can be viewed as an explanation for the way regular or chaotic behaviors 
may appear in natural systems. 

ii) Being a numerical algorithm, AS has limitations. For example, for relatively large switches 
of p or m, or for a too-large number N, A* could present some "corners". Also, obviously, 
the h size may influence the AS algorithm performances (ideally, h should decrease to zero). 
Some details and other related aspects about the errors can be found in IDanca et al., 2003$ 
and IDanca, 2003j ). 

Hi) In the general case of a dynamical system modeled by M) ; the only restriction to synthesize 
a chaotic attractor, when starting from regular attractors, is that inside the set An there are 
chaotic attractors (and vice-versa for regular synthesized attractors). 

iv) The AS can be used as a kind of control-like method IDanca, 2009$ or anticontrol f Danca et al., 20 08]. 

v) Near several continuous dynamical systems (such as the Chen system, Rossler system, 
Rabinovich-Fabrikant system (\Luo et al., 200^ ) minimal networks, Lotka-Volterra system, 
Lii system, Rikitake system), the AS algorithm was also applied successfully to systems of 
fractional orders Danca & Kai, 2010$ . 



3 Lii attractor synthesis 

The numerical results in this section are obtained using the standard Runge-Kutta algorithm 
with fixed integration time step h = 0.001. 

To visualize how the AS works, the bifurcation diagram was plotted (Fig. 3). Next, 
we synthesize the Lii attractor starting from different values for p and using deterministic 
or random schemes In this simulation, once we fixed N, all we need is to choose m 
and Vn so that the equation ^ with p* = 0.8 corresponding to the Lii attractor, can be 
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verified. Besides Poincare sections and histograms, Haussdorf distance between A* and A p * 
( |Falconer, 19 90 p. 114) was computed, in this case, in order of 1CT 2 -j- 10 -3 , which indicates 



a good approximation. 

First we applied the deterministic scheme Q for N = 2, p\ = 0.2 (corresponding to 
generalized the Lorenz system, Fig. 4 a), and p 2 = 1 (corresponding to the Chen system, 
Fig. 4 b) with the scheme [lpi,3p 2 ]- In this case, the synthesized attractor A* is identical to 
A p * with p* = 0.8 = (lp\ + 3^2) /4 (Fig. 4 c). In Fig. 4 d and e, the histograms and Poincare 
sections of both attractors, A* and A p *, are plotted superimposed to underline the identity. 

Because the solution of ^ for given N, Pn is not unique, the Lii attractor can be obtained 
in, theoretically, infinitely many ways. Thus, we chose N = 5, p\ = 0.47, P2 = 0.585, p% = 
0,678 (corresponding to the Lorenz system), p 4 = 0.905, p$ = 0.9405 (corresponding to the 
Chen system) and m x = m 2 = m 3 = m 4 = 1, m 5 = 4, again p* = 0.8 = (pi + P2 + P3 + P4 + 
4pa)/8. A Pl 5 and A*, A p * are presented in Fig. 5 with Poincare sections and histograms. 

Using the random way presented in Fig. 2, the Lii attractor can be synthesized with, for 
example, p\ = 0.6 and p 2 = 1 (Fig. 6). 

4 Conclusion 

The design AS algorithm has been utilized to generate numerically the Lii attractor starting 
from his "neighbors", the Lorenz and Chen attractors, not by continuous transformations as 
before but by discontinuous parameter switching inside a chosen parameter set. 
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Figure 1: Synthesis of a stable limit cycle for the Chen attractor, obtained using the scheme 
[lpi, 2p<2\ with pi = 0.8 and P2 = 0.959 and p* = 0.906; a) Lii attractor; b) Chen attractor; c) 
A* and A p * plotted superimposed; d) Histograms of A* and A p * superimposed; e) Poincare 
superimposed sections with plane £3 = 28 of A* and A p *. f) Time series with transients of 
component x\ of A* and A p * superimposed. 
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repeat 

label = rand(iV) 
if label = 1 then 

integrate (2) with p = p± 



inc\m 



if label = 2 then 

integrate (2) with p = p 2 



inc(m 



if label = N then 

integrate (2) with p = Pn 
inc(m' N ) 
t = t + h 
until t > T mav 



Figure 2: Random SA algorithm. 
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Figure 3: Bifurcation diagram for the Lii system. 



12 





Figure 4: The synthesized Lii attractor obtained with scheme [lpi,3p2] for p± = 0.2, P2 
1; a) A pi ; b) A p2 ; c) A* and A p * plotted superimposed; d) Superimposed histograms; 
Superimposed Poincare sections. 
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Figure 5: The synthesized Lii attractor obtained with scheme [lpi, lp2, lps, lp4,4ps] for p\ = 
0.47, p 2 = 0.585, p 3 = 0, 678, p A = 0.905 and p 5 = 0.9405. a-e) Attractors A p ., i = 1, . . . , 5; 
f) A* and A p * plotted superimposed; g) Superimpose histograms; h) Superimposed Poincare 
sections. 
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Figure 6: The synthesized Lii attractor obtained with the random scheme in Fig. 1 with uni- 
form distribution of values pi = 0.6 andp2 = La) Superimposed histograms; b) Superimposed 
Poincare sections. 
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